DMRG Study of the Ground State at Higher Landau Levels 
— Stripes, Bubbles and the Wigner Crystal — 
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Hartree-Fock theory predicted stripe or bubble phase in the third and higher Landau levels 
for two-dimensional electrons, and experimental evidences has been accumulated. In this paper 
theoretical confirmation of the stripe phase and bubble phase in higher Landau level is given by 
means of the density matrix renormalization group (DMRG) method, which can give essentially exact 
ground state for electron systems with up to 18 electrons. From the study of the pair correlation 
function, the stripe phase, bubble phase, and the Wigner crystal phase are identified, and phase 
diagram is obtained. The reentrant integer quantum Hall state is identified as the bubble state. The 
phase diagram of the fourth Landau level shows more diversity than the third level. 
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I. INTRODUCTION 



In the lowest Landau level the fractional quantum Hall 
effect is observed. However, it has not been observed in 
the Landau levels with Landau index ATl larger than two. 
This is because the interaction between electrons is not 
short-ranged enough to make Laughlin state realized. In 
such higher Landau levels, Koulakov et al. proposed re- 
alization of the stripe and bubble phases based on the 
Hartree-Fock theory [1,2]. These phases belong to fam- 
ilies of the charge density wave states proposed before 
the discovery of the fractional quantum Hall effect in the 
lowest Landau level [3,4]. So it is natural to anticipate 
such states once the FQH states are denied. 

The evidence of the stripe phase has been obtained 
from experiments on ultra high mobility samples at low 
temperature as remarkable anisotropy in the longitudi- 
nal resistance when the higher Landau level is nearly 
half-filled [5,6]. Although no direct evidence of the bub- 
ble phase has been obtained, the reentrant IQHE around 
v* ~ 1/4 and 3/4 [7] suggests the realization of the bub- 
ble phase, where v* = v — [v] is the filling factor of the 
topmost partially filled Landau level. 

Since, the Hartree-Fock theory may fail to find the true 
ground state as evidenced in the lowest Landau level, 
more reliable theory is desirable. So exact diagonaliza- 
tion studies have been done, and the realization of the 
stripe phase and the bubble phase in the ground state 
was confirmed [8]. However, the system size treated by 
the exact diagonalization is quite small, and the filling 
factor studied was limited. Because of these limitations, 
such study cannot determine the phase diagram of the 
present system. In the present study we investigated the 
ground state by DMRG method [9]. We have investi- 
gated systems with up to 18 electrons, and determined 
the ground state at various filling factors by examining 
the pair correlation function. Part of the investigation, 
where we clarified the phase diagram in the third lowest 



Landau level (N^ = 2) has already been published [10]. 
Here, we give new results for the fourth lowest landau 
level together with the previous results. 



II. MODEL AND METHOD 

The Hamiltonian for electrons in Landau levels con- 
tains only the Coulomb interactions. After the projection 
onto the A^th Landau level, the Hamiltonian is written 
as 



(1) 



where is the guiding center coordinate of the zth elec- 
tron, which satisfy the commutation relation [i?J , i?^] = 
iSj t k- Ln 1j (x) are the Laguerre polynomials, and V(q) = 
2ne /q is the Fourier transform of the Coulomb interac- 
tion. The magnetic length I is set to be 1. We consider 
uniform positive background charge to cancel the com- 
ponent at q = in eq. (Q), and neglect the electrons in 
fully occupied lower Landau levels. Periodic boundary 
condition is imposed both in the x and y directions the 
period being L x and L y , respectively. The filling factor 
of the relevant Landau level is v* = 2nN c / L x L yi where 
N c is the total number of electrons. We used the DMRG 
method to obtain the ground state of the system with 
up to 18 electrons at various fillings with changing the 
aspect ratio L x /L y . 

We first investigated the aspect ration dependence of 
the ground state energy at fixed filling factors. Then pair 
correlation function g(r) in guiding center coordinates, 



L X Ly 



is calculated at local minima of the ground state energy. 
For the most filling factors our system is large enough 
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that it is rather easy to tell the ground state CDW pat- 
tern by looking at g(r). 



III. RESULTS 

From the pair correlation function we found that the 
ground state at Ni, = 2 and 3 is either the Wigner crys- 
tal, the bubble state, or the stripe state depending on the 
filling factor v* and the Landau level index TVl, except for 
a small region of the filling factor around 0.28 in N-^ = 3 
Landau level. The resultant phase diagram is shown in 
Fig.l and Fig. 2. For Nl = 2 Landau level there are three 
phases. The boundary between the Wigner crystal state 
and the two-electron bubble phase lies between v* = 2/9 
and 1/4, and that between the two-electron bubble and 
the stripe phase lies between v* = 4/11 and 2/5. In- 
spection of the r dependence of g(r) shows that these 
boundaries are the first order phase transition points. 
For Ml — 3 Landau level there are at least four phases. 
Here three-electron bubble phase, where cluster of tree 
electrons form a kind of the Wigner crystal, becomes 
possible. The ground state between v* — 1/4 and 0.3 
has not been identified. We will discuss this state later. 
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FIG. 1. The phase diagram for the third lowest Landau 
level (7V L = 2) obtained by DMRG and Hartree-Fock meth- 
ods. 
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FIG. 2. The phase diagram for the third lowest Landau 
level (Nh ~ 3) obtained by DMRG and Hartree-Fock meth- 
ods. 



We also solved the Hamiltonian by the Hartree-Fock 
theory [4]. In this method the filling factor is also lim- 
ited to simple rational fractions. However, for the bub- 
ble phase the order parameters obtained by the this 
method almost coincide with those proposed by Fogler 
and Koulakov [2], just the same as the case of the Wigner 
crystal state [4] . Then it is possible to obtain the Hartree- 
Fock ground state as a continuous function of the fill- 
ing factor, and transition point between different bubble 
states can be determined. The stripe phase has weak 
modulation along the stripes, which opens the gap at the 
Fermi level and lowers the energy. The phase diagram by 
this method is also shown in Figs.l and 2. The phases 
are separated by the first order phase transition similarly 
to the DMRG case. The phase diagrams from these two 
methods are similar except that the correct phase dia- 
gram by the DMRG shows wider region for the stripe 
phase. The three-electron bubble phase at A^l = 2 and 
the four electron bubble phase at Al = 3 are swallowed 
by the stripe phase. 



IV. DISCUSSION 

In this paper we have determined the phase diagram 
of electrons in the third and fourth lowest Landau lev- 
els by the DMRG method. The slight difference in the 
phase diagrams for A^l = 2 and 3 comes from the range of 
the Coulomb interaction projected to each Landau level, 
namely the Laguerre polynomial in eq.(|l|) is responsible 
for the difference. In actual systems there are several 
factors which affect the form of the Coulomb potential. 
The first is the spread of the wave function in the direc- 
tion perpendicular to the two-dimensional plane. This 
spread acts mostly to reduce the short-range part of the 
Coulomb interaction. Since this part is already reduced 
by the Laguerre factor in the higher Landau levels, we do 
not expect this will cause qualitative change to our phase 
diagram. Another factor is the screening by the electrons 
in the filled lower Landau levels. This acts to reduce the 
long-range part of the interaction. We did not consid- 
ered this factor, since this effect is sample dependent. 
Namely, it depends on the electron density as well as the 
filling factor. Comparing our Hartree-Fock results with 
those by Shklovskii's group, where the screening effect 
is considered [1,2], we have concluded that the screening 
effect does not affect the qualitative feature of the phase 
diagram. 

We think that the reentrant IQH state is the two- 
electron bubble phase, which is pinned by impurities. 
Then our calculation predicts that the reentrant phase 
should be observed in ATl = 3 Landau level at lower fill- 
ing factor. 

Our calculation until now have not identified the 
ground state between v* — 1/4 and 0.3 for A/l = 3 Lan- 
dau level. The g(r) looks like stripe phase with stronger 
modulation. Actually the energy of stripe phase by the 
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Hartrec-Fock theory is quite close to those of the 2- 
electron and 3-electron bubble phases here. Since the 
energy lowering by the quantum fluctuation is larger in 
stripe phase as the result around v* ~ 1/2 shows, reen- 
trant stripe phase is a feasible candidate. We hope that 
we can identify this phase in near future. 



ACKNOWLEDGMENTS 



Part of the numerical calculation is performed in the 
ISSP, University of Tokyo. D.Y. thanks Aspen Center for 
Physics, where part of the work was done. The present 
work is supported by Grant-in-Aid No. 12640308 and 
No. 11740184 from MEXT, Japan. 



[1] A. A. Koulakov, M.M. Fogler, and B.I. Shklovskii, Phys. 



Rev. Lett. 76, 499 (1996), and M.M. Fogler, A. A. 
Koulakov, B.I. Shklovskii, Phys. Rev. B 54, 1853 (1996). 
[2] M.M. Fogler and A. A. Koulakov, Phys. Rev. B 55, 9326 
(1997). 

[3] H. Fukuyama, P.M. Platzman, and P.W. Anderson, 

Phys. Rev. B 19, 5211 (1979). 
[4] D. Yoshioka and H. Fukuyama, J. Phys. Soc. Jpn. 47, 

394 (1979), D. Yoshioka and PA. Lee, Phys. Rev. B27, 

4986 (1983). 

[5] M.P. Lilly, K.B. Cooper, J.P Eisenstein, L.N. Pfeiffer, 
and K.W. West, Phys. Rev. Lett. 82, 394 (1999). 

[6] R.R. Du, D.C. Tsui, H.L. Stormer, L.N. Pfeiffer, K. W. 
Baldwin, and K. W. West, Solid State Commun. 109, 
389 (1999). 

[7] K.B. Cooper, M.P. Lilly, J.P. Eisenstein, L.N. Pfeiffer, 
and K.W. West, Phys. Rev. B 60, R11285 (1999). 

[8] E.H. Rezayi, F.D.M. Haldane, and K. Yang, Phys. Rev. 
Lett. 83, 1219 (1999), F.D.M. Haldane, E.H. Rezayi, and 
K. Yang, Phys. Rev. Lett. 85, 5396 (2000). 

[9] S.R. White, Phys. Rev. Lett. 69 (1992) 2863; Phys. Rev. 
B 48 (1993) 10345. 
[10] N. Shibata and D. Yoshioka, Phys. Rev. Lett. 86 No.25 
(2001). 



3 



